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^ ; Abstract 

We study a softly-broken supersymmetric model whose gauge symmetry is that of the standard 
JL , model (SM) gauge group times an extra Abelian symmetry U(lY. We call this gauge-extended 

model U(l)' model, and we study a U(l)' model with a secluded sector such that neutrinos acquire 
Dirac masses via higher-dimensional terms allowed by the U(l)' invariance. In this model the jit term 
of the minimal supersymmetric model (MSSM) is dynamically induced by the vacuum expectation 
value of a singlet scalar. In addition, the model contains exotic particles necessary for anomaly 
cancellation, and extra singlet bosons for achieving correct Z'/Z mass hierarchy. The neutrinos are 
charged under U(l)' , and thus, their production and decay channels differ from those in the MSSM 
in strength and topology. We implement the model into standard packages and perform a detailed 
analysis of sneutrino production and decay at the Large Hadron Collider, for various mass scenarios, 
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concentrating on three types of signals: (1) 0£ + MET, (2) 21 + MET, and (3) U + MET. We 
compare the results with those of the MSSM whenever possible, and analyze the SM background 
for each signal. The sneutrino production and decays provide clear signatures enabling distinction 
of the U(l)' model from the MSSM at the LHC. 

PACS numbers: 12.60.Cn,12.60.Jv,14.80.Ly 
Keywords: Supersymmetry, Scalar Neutrinos, LHC 
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I. INTRODUCTION AND MOTIVATION 

The MSSM is arguably the most popular 'new physics' scenario referring to a perturbative 
completion of the SM beyond Fermi energies. Motivated by the resolution of such long 
standing problems of the SM as the gauge hierarchy problem, the existence of dark matter 
and the added attraction of gauge unification, nevertheless, it still has some outstanding 
problems. One of these is the so-called fj, problem (l| . Supersymmetric models which extend 
the MSSM via an extra gauge group generally intend to solve // problem and incorporate an 
extra singlet field, whose coupling to the Higgs fields and VEV generate dynamically the \i 
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term. These models extend the SU(2)l <g> U(1)y MSSM electroweak symmetry by an extra 
U(l) gauge symmetry. Such an extension is minimal, and it is well motivated in superstring 
theories 2|, grand unified theories [3] and in dynamical electroweak breaking theories J4|. 
The simplest versions contain a singlet field and an extra neutral gauge boson. Other 
versions also allow right-handed neutrinos into the spectrum. In a non-minimal version of 
the U(l) extended MSSM, which includes several singlet (S) fields, the tension between the 
electroweak scale and developing a large enough Z' mass is resolved. We call this version of 
the model secluded sector U(l)', a shorthand notation for SU(3) c ®SU(2)l®U{1)y®U(1)', 
the gauge symmetry underlying the model, and describe it in the next section. In the MSSM, 
as in the SM, neutrinos are massless. The fact that neutrino oscillation imply non-vanishing 
neutrino masses is a strong motivation to consider an extended form of the MSSM. Small 
neutrino masses consistent with neutrino oscillation phenomenology are usually explained 
by the see-saw mechanism |5|. In the see-saw mechanism, large Majorana masses for right- 
handed neutrinos induce small Majorana masses for left-handed neutrinos. In the scalar 
sector, right-handed sneutrinos mix with the left handed sneutrinos and give potentially 
new signals for extended symmetry. The choice of U(l) symmetry would determine the 
magnitude and type of neutrino masses. In this paper, we consider a U(l)' extended form 
of the MSSM that contains Dirac-type neutrino masses. 

Direct or indirect detection of the superpartners of the Standard Model particles, the 
definitive signal for supersymmetry, is one of the major aims of the LHC experiments. Except 
for the LSP in the R— parity conserving supersymmetry, the superpartners are expected to 
decay instantaneously into SM particles, plus the LSP, detected as missing energy. The 
common methodology for detection is to analyze the production and cascade decays of the 
supersymmetric particles. As the right sneutrinos, which can mix with the left sneutrinos 
are a feature of the U(l)' model that distinguishes it from MSSM, studying sneutrino signals 
would be an important test for this model. 

Systematic analyses of sneutrino decays in the MSSM have been performed in jf|. The 
aim of this article is to perform a comparative study of LHC signals of sneutrino production 
and decays in the MSSM and in a supersymmetric model with a secluded U(l)' breaking 
sector Q via their decay chain topologies. Differences between MSSM and the secluded 
sector U(l)' model likely reveal themselves via decay modes of the sneutrino. We analyze 
the signals, and, for completeness, we also include possible Standard Model backgrounds. 



In most variants of the MSSM consistent with relic density calculations, the LSP is 
the lightest neutralino, typically a mixed state of bino (fermionic partner of the U(l)y 
gauge boson) and the higgsino. In a previous work [8], we showed that a minimal U(l)' 
model (one extra singlet boson) could be consistent with the excess positron observed in 
satellite experiments, choosing on of the right-handed sneutrinos as the LSP. However, for 
the purpose of this work (dependent on parameter space chosen to compare our results with 
those of MSSM), the secluded sector U(l)' lightest neutralino appears consistently to be 
the lightest supersymmetric particle (LSP) and therefore is a potentially viable dark matter 
(DM) candidate, although its composition is likely to differ from the lightest neutralino in 
MSSM. 

Here we perform a thorough analysis of sneutrino production and decay in the secluded 
sector U(l)' model. In order to compare with previous signals, we establish a set of three 
mSUGRA-inspired benchmarks for our model. Similar to the mSUGRA benchmark points 
analyzed in MSSM (LM1, LM2, LM6) 9 n||, we analyze the corresponding scenarios in 
secluded sector U(l)' model (LM1', LM2', LM6'). Here LM stands for Low Mass, a choice 
likely to yield visible signals at the LHC 

Our paper is organized as follows. We briefly introduce the model in Section [Til then define 
the parameters and physical masses of supersymmetric particles in the secluded sector 11(1)' 
model in Section [TTTI For each benchmark point, we insure that DM candidate of the U(l)' 
model yields relic densities consistent with the WMAP range of cold dark matter density 



12fl . We then perform a comparative analysis of the production, decays and detectability 



of sneutrinos within these benchmark supersymmetric scenarios. During this analysis we 
focus on the multilepton plus missing energy signatures of the supersymmetric scenarios. 
We present the results of our simulation analysis for the LHC. In Section IIVI we conclude 
the work. We leave the extensive details of the model for the Appendices. 



II. THE U(iy MODEL 

The MSSM suffers from a naturalness problem due to the presence of fi parameter, 
responsible for giving masses to the Higgs bosons and Higgsino in the superpotential. From 
a purely theoretical point of view, the value of this parameter is expected to be either of the 
order of the GUT, Planck scale or zero. For phenomenological aspects, however, it must be 
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of the order of the scale of electroweak symmetry breaking (EWSB) and it has to be non-zero 
to agree with the experimental data. Seen from the low energy point of view, adding an 
extra U(l) is needed in order to solve the // problem {l| of the MSSM. Basically the problem 
is remedied by extending the matter and gauge structure of the MSSM, e.g. within unified 
and/or string models by introducing an additional singlet filed S, whose VEV generates the 
fi term dynamically. Theories with an extra U(l)' broken at the electroweak- to- TeV scale 
by SM singlets are known to be able to generate an appropriately sized \i parameter (see 

e.g. fl). 

The other success of the U(l)' symmetry is being able to generate pertinent neutrino 
masses by introducing right-handed neutrinos into the superpotential. The right-handed 
neutrino sector and the /i parameter can be correlated for both Majorana [13| and Dirac 



masses 14]. We assume here that lepton number is an accidental symmetry that is conserved 
at the perturbative level. Hence, the neutrinos are Dirac fermions, requiring Yukawa cou- 
plings of O (10 -13 ). These couplings are technically natural, but an explanation for such a 
strong suppression is clearly desirable. One way this can occur is if the U(l)' invariance sup- 
presses leading order contributions to Dirac neutrino masses and allows higher- dimensional 

n 

operators [14]. 

In this work, we extend the MSSM in the following ways. First, the gauge structure of the 
MSSM, SU (3) SU (2) l®U (l)y , is enriched to include an extra Abelian group factor U(l)'. 
Second, we promote the fi parameter into dynamical field, S, which is charged under the 
U(l)'. Third, exotics with Yukawa couplings to S are included to make the theory anomaly- 
free. Fourth, Z'/Z mass hierarchy in the model is ensured by three additional SU(2) singlet 
fields which are coming from secluded sector of the model. The model also includes a term 
that provides suppressed Dirac neutrino masses in accordance with observations. We present 
the main relevant points in this section, leaving the details for the appendices. 

In the minimal version of the model which contains only one singlet S, there is some 
tension between the electroweak scale and the need to generate a large enough M' z . These 
two problems can be decoupled without fine tuning when several additional fields are incor- 
porated into the model. An example of this kind of non-minimal model is secluded sector 
model. The secluded sector model involves an ordinary sector of symmetry breaking fields, 
which includes two Higgs doublets, and an SU(2)l singlet S. After acquiring a VEV, S 
generates an effective /j parameter /j = h s (S). The secluded sector of the model includes 



Field Q U D L N E H u H d S S\ S2 S~ 3 Q Q £ £ 

SU(3) C 3331111 111113311 

SU{2) L 2112112 211111 111 

U(1) Y 1/6 -2/3 1/3 -1/2 1 1/2 -1/2 Y Q -Yq Y c -Y c 

^(i) Qq Q'u Q'd Q'l Q'n Q'e Q'h u Qiid Q's Qsi Qs 2 Qs 3 Q'q Q'q Q'c Q'z 

TABLE I. Gauge quantum numbers of quark (Q,U,D), lepton (L,N,E), Higgs (H u ,H d ), SM- 
singlet (S, S\, S2, S3), exotic quark (Q, Q) and exotic lepton (£,£) superfields. 

three SU(2)l singlet fields S^i = 1,2,3 which acquire large VEVs. All four VEVs of the 
singlet fields S,Sx t 2,3 contribute to Z' mass. Thus, in this model, Z'/Z mass hierarchy is 
implemented mainly through the secluded sector of the model. 
The superpotential of the model is given by 

W = h u Q- H U U + h d Q ■ H d D + h e L ■ H d E + h s SH u ■ H d + -^—S\L ■ H u h u N + hS^Ss 

Mr 

+ ^h^SQM + y £hisc j % (1) 

i=l j=l 

where the fields entering the equation, together with their quantum numbers are listed in 
Table HI Here, Mr is a large mass scale and h v is the Yukawa coupling responsible for 
generating neutrino masses. 

The U(l)' charges of the fields satisfy a number of conditions arising from phenomeno- 
logical constraints, as well as from gauge invariance of the model and from the requirement 
of cancellation of gauge and gravitational anomalies. They are as follows. 

The U(l)' charges satisfy Q'h u + Q'u A 7^ to forbid the bare fi term, Q' L + Q'h u + Q'n 7^ 
to induce neutrino masses correctly, and Q' Si + Q'g 2 + Q' Sa = to correctly generate the 
Z — Z' mass hierarchy. Gauge invariance of the superpotential implies 

= Q's + Q' Hu + Q' Hd , 
o = Q' Q + Q' Hu + Q'u, 
o = Q' Q + Q' Hd + Q'D, 

= Q'l + Q'n d + Q'e, 
= Q' q + Q'q + Q' s , 
= Q' c + Q' z + Q's, 
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o = Q , Sl + Q , L + Qk + Q , N- (2) 

For the model to be anomaly-free the U(l)' charges of fields must satisfy 

= 3(2Q' Q + Q' v + Q' D ) + n Q (Q' Q + Q'^) , (3) 
= 3(3Q' Q + Q' L ) + Q' Hd + Q' Hu , (4) 

o = 3(^g' Q + \q' d + iqf v + ^g' L + q' e ) + ~(g^ + g' H j 

+ 3n Q Y 2 (Q' Q + gy + n c Y 2 (Q' c + QL) , (5) 
= 3(6g' + 3Q' V + 3Q' D + 2Q' L + Q' E + Q' N ) + 2g'^ + 2Q' Hu 

+ Q's + Qs a + Q's 2 + Q's 3 + 3n s (Q' Q + gy + n £ (g^ + Q' z ), (6) 
= 3(g' Q 2 + Q'l - 2Q' 2 - Q' L 2 + Q'i) - Q'l + Q' H \ + 3n Q Y Q (Q' Q 2 - QL 2 ) 
+ n c Yc(Q' c 2 -Q' z 2 ), (7) 



o = 3(6g' 3 + 3g' D 3 + sg'^ 3 + 2g' L 3 + g' E 3 + g^ 3 ) + 2g' H 3 + 2g' H 3 + g' 5 3 

/ 3 , r\> 3 



+ Q's! + Q's! + Q's! + 3^ Q (g' Q 3 + Q'q 3 ) + n c (Q' c 3 + QL?), (8) 

which correspond to vanishing of U(l)'-SU(3) c -SU(3)c, U(1)'-SU(2) L -SU(2) L , U(l)'- 
U(1)y-U(1) y , [/(l)'-graviton-graviton, U(1)'-U(1)'-U(1) Y , and U {!)' -U {!)' -U {!)' anomalies, 
respectively. All these anomaly cancellation conditions are satisfied for a particular pattern 
of charges and parameters. It is found that the solution to the mixed anomaly constraints 
requires hq = 3 color triplet pairs with hypercharge Yq = —1/3, and nc = 5 singlet pairs 



with Yc = —a/2/5. With these parameter values one obtains the U(l)' model displayed in 
Table HT1 The U(l)' charges for Higgs fields in the model are chosen as 

Q's = -Q' Sl = -Q's, = \Q's 3 , Q' Hu + Q' Hd + Q' s = o. (9) 

Under the conditions above, the supersymmetry breaking soft terms for the secluded sector 
model are 

Vsoft = Ko/t + v s ° oft (io) 

where Vj a f t are the allowed U{1)' dimension-2 operators 

VLft = (m 2 SSl SSt + m 2 SS2 SS 2 + m 2 SlS2 S{S 2 + h.c.) (11) 

and V° a f t term is defined as 

Kft = m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 \S\ 2 + (12) 
- (A s h s SH u H d + A,h s S l S 2 S 3 + h.c.) 



Q'h u = "2 

Q' Hd = i 

Q' s = l 

Q' Sl = -i 



Qk = - 

Q's 3 = 2 
Q'q = x 

Q'u = 2- 



Q'd 
Q'l 



— 3x 



Q E — 3 



Q' Q = 

Q'c = 



18 

-22+12^+^^ 
18 

30 

-15-13v / lQ+12y / l0a:-v / 5n 
30 



TABLE II. A set of U (1)' charges satisfying all gauge invariance and anomaly cancellation condi- 
tions. The charge of the quark doublet Q is left free, and for simplicity £l(x) = a/241 + 708x + 612x 2 
is introduced. 

We set ttl 2 SiS2 = as only two of the Si fields are needed to break the global U(l) symmetries. 
To insure that the potential is not bounded from below, we require 

m 2 s + m 2 Sl + 2m| Sl > , 
m 2 s + m 2 S2 + 2m 2 s S2 > °- 



(13) 



In the model, the charge of the quark doublet Q is kept as a free parameter after the 
normalization Q' Hu = -2, Q' Hd = 1, Q' s = 1, Q' Si = -1, Q' S2 = -1, Q' Sa = 2. 

In this model the left and right sneutrinos mix, and the mixing matrix can in general be 
expressed as 



(ml m 2 tJ 



LB 



1 1 

mtij mtj 




(14) 



where i, j are the flavor indices and the matrix elements are given by 



ml lL = Ml + (m-) 2 + - A (g 2 Y Y L - »-)((H») 2 - (H°) 2 ) 
1 



+^y'Q'l(q'hM 2 + q'hMY + Q's(s) 2 p s ) 



ml, = M 2 Nj + «)' + -g Y Y N ((H^ - (H") 2 ) 

RR J 4 

+ 1 t;9y'Q'n(Q / h u (H u } 2 + Q' Hd (H° d } 2 + Q's(S) 2 Ps) 



2 

( m h ) 



m 



a* P , h s (S 2 ){S 3 ) 

A »i + + 



(15) 



tan/3 V2(S!) 

Here M|. and Mjy. are soft mass terms and A v . are trilinear couplings (assumed flavor- 
diagonal). Dirac neutrino masses m u , the fi parameter and p s in the equations above are 
expressed as 

Ps = i + E h Q 'f l - as) 



M, 



R 



V2 ' 



Q's v s 



In this model, neutrino masses are chosen to be Dirac-type. The effective neutrino Yukawa 
coupling Y„ leads to neutrino masses in agreement with the experiment. Numerically, we 
obtain [8] 

12 \ 1/2 



3 x 10 



-13 



2.8 x 10" 3 eV 



(17) 



III. MSSM VS. £/(!)' AT THE LHC ENERGIES 



A. Parameter Space and Relic Density 



Motivated by the fact that the scalar neutrino LSP can in principle explain the WMAP 
data as well as excess positron flux measured by various satellite experiments |8|], we ana- 
lyze the model further by investigating the production and decay mechanism of the scalar 
neutrinos at LHC. 

The model we consider here, the sechided-{7(l)' with right-handed neutrinos, has some 
advantages over the so-called the minimal U(l)' where only one additional scalar field is 
introduced. The squark phenomenology in this minimal U(l)' model has been explored in 
Ref. (ij where there is difficulty with inducing a small /i e // while satisfying the Z' mass 
bound, which is around 1 TeV. This is because both fi e ff and mz> are proportional to 
the vacuum expectation value of the additional scalar field S. One needs three additional 
scalars to ameliorate the picture the VEVs of the new scalars are kept large large. This is 
one motivation for the secluded U(l)' model. For further details of the model, see 




V L.R 



"L.R 




H\ . H2 , . . . , H 6 



V L.R 



V L.R 



FIG. 1. The Feynman diagrams for the production of the scalar neutrinos in the secluded-U {!)' 
model. Hi,i = 1, ...,6 are the CP-even physical Higgs bosons. 



The Feynman diagrams contributing to the hard production of scalar neutrinos are given 
in Fig[TJ For simplicity we neglect the mixing between Z ( the Z boson of the SM) and Z' in 
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the numerical analysis so that Z\ = Z and Z2 = Z'. In addition to Z' exchange (left-handed 
diagram), all CP-even Higgs bosons contribute to the process in the s-channel (right-handed 
diagram) . 

Once the scalar neutrinos are produced, they will decay. The decay pattern strictly 
depends on scenario chosen for the free parameters. Since we are interested in rather light 
scalar neutrinos (assuming low-energy SUSY exits), we prefer to choose MSSM-like Low 
Mass (LM) scenarios Battaglia et. al have proposed updated post-WMAP benchmark 
points for the constrained MSSM [10] modifying earlier proposal fll| . and we include these 
points in Table IIHI 

To compare our results with MSSM predictions, we choose three low-mass MSSM sce- 
narios (benchmark points), namely LM1, LM2 and LM6, from the low mass scenarios of 
MSUGRA and use Softsusy package [16| to generate the MSSM spectrum. In the se- 
cluded U(l)' we choose LM-like scenarios, denoted as LM1', LM2' and LM6' by keeping 
the overlapping parameters the same and fixing the additional parameters to agree with 
phenomenological constraints on masses. The input parameters for LM1, LM2 and LM6 
for MSSM as well as their corresponding prime versions for the secluded U(l)' are given 
in Table IIHI As seen from Table IHIt the VEVs of the additional scalars (Si,S% and S3) 
v Si ,i = 1, 2, 3 are taken above the TeV scale so that the Z' mass bound is satisfied no matter 
what the VEV of the scalar field S is chosen. In fact, for convenience, the parameters /j, e ff 
and h s are taken as free parameters and the VEV of S are determined accordingly using 
the relation given in Eq. (TT6j) . From Table [Tvl it is seen that the scalar neutrino masses are 
rather light. The left-handed sneutrinos masses are varying in the 168 GeV-287 GeV range 
while the right-handed ones are in the 412 GeV-704 GeV depending on the LM scenario as 
well as on the flavor of the scalar neutrino. The right-handed scalar neutrinos are heavier, 
showing the same pattern as in the neutrino sector. With these chosen masses we can foresee 
that the production cross section for the left-handed sneutrinos will dominate the one for 
the right-handed ones. 

The validity of the MSSM scenarios LM1, LM2 and LM6 has been confronted with both 
the LEP and Tevatron data. There will be no contributions to the LEP observables from 
our LMX', X = 1, 2, 6 scenarios since the lightest Higgs boson mass in the model is 218 
GeV, which is already above the LEP energy. For the Tevatron case, however, one needs 
do a more careful analysis. Nevertheless, as the LMX' scenarios aim to be consistent with 
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TABLE III. The scenarios (benchmark points) LM1, LM2, and LM6 (for the MSSM i. e. minimal 
super gravity), and LM1' , LM2' and LM& (for the U(l)' model). The unprimed LMX and primed 
LMX' benchmark points similar mass spectra. Parameter Ry is defined in AppendixlCH 
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TABLE IV. The complete mass spectra of the benchmark points (scenarios) given in Table [7711 for 
both MSSM and the secluded U{1)'. 



the corresponding MSSM scenarios, in the limit where the extra U(l)' particles decouple, 
we expect consistency with the Tevatron data. To verify this point, we used the package 
HiggsBounds 1^|, which yields results for any arbitrary Higgs sector. 
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TABLE V. The production cross section and the relic density Odm values for the LM scenarios 
considered in the paper. 



The production cross sections for the scattering pp — > ^e L R ^e L R processes are listed 
in Table [V] for both MSSM and the secluded U(l)' model. The values were obtained 
implementing the secluded U(l)' model into CalcHEP [l8| with the help of LanHEP [l9|. The 
parton distributions in the proton have been parametrized by using CTEQ6M of LHAPDF ^20 1 . 
The MSSM total cross sections (including the three scalar neutrino flavors) are in the range 
of 4 to 110 fb while in the secluded U(l)' model they are varying between 1.1 pb to 2.6 
pb. The new right-handed sneutrino cross sections in the secluded U(l)' model are about 
10 times smaller than the cross sections for their left-handed counterparts, and are in the 
range of 20 fb to 80 fb. 

In Table |V] we also included the relic density of the dark matter for all six scenarios. 



This calculation is straightforward using the Mi cr omegas package 



2l| . once we include the 



model files from CalcHEP. All the numbers are within the la range of the WMAP result 



12] 
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which can be given with those from the Sloan Digital Sky Survey 

Q DM h 2 = 0.1111°;°^. (18) 
We note that the relic density of the dark matter QdmIi 2 is very sensitive to the free parameter 
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FIG. 2. The Feynman diagrams for the three-body decay channels of the next-to-LSP X2 Jri the 
secluded-U (1)' model. Here Hi,i = 1,...,6 are the CP-even physical Higgs bosons while A®,i = 
1, ...,4 are CP-odd ones. 

Ry listed in Table HTT1 which varies between 42 to 50. It's defined (see also Appendix [C]) as 
the ratio between bare U(l) gaugino masses 

R Y > = My, I My 

where M Y and My, are the Bino and Bino' mass parameters appearing in the 9x9 neu- 
tralino mixing matrix. More details are given in Appendix In Table [IV] the Lightest 
Supersymmetric Particle (LSP) is the lightest neutralino Xi with masses 96 GeV, 63 GeV 
and 79 GeV for the LM1', LM2' and LM6' scenarios, respectively. The next-to-lightest su- 
persymmetric particle is X2 with masses 99 GeV, 138 GeV and 158 GeV, respectively. For 
such a spectrum, there will be no kinematically available two-body decays for the xHi so 
that three-body channels need to be considered. The three-body decay modes relevant to 
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the analysis here are given in Fig. [2j These decay modes will be considered in the LHC 
simulation but not the relic density calculation, as they give negligible contributions. The 
sizable contributions to the relic density are for the LM1' scenario 

• XiX^t-t+ (15%) 

• Xlxl^r-r+ (13%) 

• Xi X2 ~ e~e + / / u~/i + (8% for each channel) 

• X2X2 ~^ e~e + /{i~fi + (8% for each channel) 

• X2 X2 ~ * v i^i ,1 = e i^) T (5% for each channel) 

• x\x°i^t-t+ (5%) 

• Xi Xi ~ > e~e + /fi~fi + ( 3% for each channel) 

• X1X2 ~^ v i^i 1 1 = e i ^1 T (3% f° r each channel) 

• X°iXl^W-W+ (2%) 

In the LM2'(LM6') only Xi Xi annihilation contributes to the relic density of the dark matter 
as follows 

• X?X?->r"r+ (75% (38%)) 

• X?X?-^> + (8% (26%)) 

• X?X?^e-e+ (8% (26%)) 

• X°iX°i^bb (3% (%1)) 

• xSxS-^frPr (1%(3%)) 

• X°iX°i^mJ = e,fi (0% (3%)) 

• Xlxl^dd/ss (1% (0%)) 
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Contributions from the X2 X2 or X? X2 annihilations for the Xi scenario are due to the fact 
that Xi an d X% are almost degenerate in mass and since the two-body decay channel limit 
is used in Mi cr omegas, X2 ac ts very similar to Xv There is no sizable contributions from X2 
in the other two scenario since X2 * s m uch heavier. The bino, wino, higgsino and singlino 
compositions of the neutralinos for the scenarios LM1', LM2' and LM6' are given in Table [VTl 
in Appendix [D] The LSP Xi is mainly bino for LM1' but mostly singlino (S) for LM2' and 
LM6' (94.2% and 93.6%, respectively, for the two scenarios). The situation is reversed for 
the next-to-LSP, Xt 

B. The LHC Signals 

After discussing the chosen scenarios and the details of the relic density calculation of 
the dark matter, we proceed to discuss the signals at LHC from scalar neutrino production 
processes. To determine and classify all possible signals for the scenarios LM1', LM2' and 
LM6', we need to look at the decay topology of the scalar neutrinos. 

Since we include MSSM scenarios LM1, LM2 and LM6 for comparison purposes, we first 
outline the main decay modes governing the decay channels. The left-handed scalar neutrinos 
i>i L decay to Vg x\ with about 100% branching ratio for the LM1 and LM6, since all the other 
neutralinos are heavier than the scalar neutrinos. The picture is a bit more complicated for 
the LM2 where v iL ,£ = e,/i decay to v iX \ (71%), Ixt (20%) and v lX \ (8.8%). For the u TL , 
the branching decay ratios are Wt\ (61.5%), v T Xi (34.2%), rxf (3%) and v T x° (1.3%). 
Further in the decay chain X2 decays mainly to tti/tt^ (48% for each channel), and the 
chargino x^tov T T\ (with 95.4% branching ratio) and W Xi (4.6% branching ratio). 

In the secluded U(l)' model, the decay modes of the scalar neutrinos with more than 1% 
branching ratio are, for the scenarios LM1'/LM2 / /LM6 / 

• %(%J -+v t )&, 8.6% (0%) / 91.7% (84.8%) / 93.2% (65.7%) 

• v ltt (v lR ) -+u t )$, 91.4% (90.6%) / 8.3% (0%) / 6.8% (0%) 

• v iL (u tR ) ^vix!, 0% (8.8%) / 0% (14.8%) / 0% (34.1%) 

There will be further decays of X2 an d Xa i n the chain. It is better to consider xi first. 
Again in the same notation (LM1'/LM2'/LM6 / ) it decays as 
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• xl -> tt* (tt*) , 6.5% (4.9%) / 9.1% (3.7%) / 8.6% (3.9%) 

• xl^^l (th) > 6 -°% ( 5 - 3 %) / 5-2% (5.0%) / 4.0% (8.3%) 

• xl^^K^KJ, 5.2% (5.2%) / 5.4% (5.9%) / 4.2% (4.3%) 

where £ — e, fi and the conjugated decay modes are not listed. Then the decay modes of the 
scalar leptons for £ = e, are 

•h~>£Xi (^X§) , 6.7% (93.3%) / 90.1% (9.9%) / 92.6% (7.4%) 

• Zr -> ^X? (€5g) , 20% (80%) / 71.6% (28.4%) / 89.7% (10.3%) 

As can be seen from the above decay patterns, each decay ends up with either an LSP 
Xi or next-to-LSP particle \\. As mentioned earlier, \\ cannot decay into two-body but 
instead must undergo the one of the three-body decays given in Fig. [2j The relative ratios 
arc| given in the (LM1'/LM2 / /LM6 / ) order as 

• X2 vt »i {y T v T ) x? , 24% (24%) / 4.5% (6.3%) / 1.8% (1.8%) 

• x°2 f + Z~ (r + T~) Xi , 14% (0%) / 10.8% (63%) / 28.7% (37%) 

In the light of these decay patterns, there are mainly three types of signal: (1) 0£+ $t, 
(2) 2£+ $t and (3) 4£+ It is i n fact also possible to produce signals with six or eight 
leptons but the probability is very suppressed thus we ignore such signals. Therefore, in the 
rest of this section we discuss these three signals at LHC. Predictions from MSSM will be 
included as well. In MSSM there is no the 4£+ type of signal in MSSM for the LM1 and 
LM2 and LM2 scenarios. The 2£+ $t signal is possible through chargino xf decay. 

The usual concern is the possible background for the signals from the SM. For the 0£+ $t 
mode, the background will come from the Drell-Yan (D-Y), pp — > U£p£, and pp — > ZZ where 
each of Z decays invisibly. Since the D-Y has a huge cross section, some cuts need to be 
implemented. In the 2£+ ft? case, in addition to the D-Y and ZZ production (where 
one of the Z decays leptonicaly) , there is W + W~ production. In principle there could be 
contributions from the tt with jet veto, but we ignore such possibility since the 6-jets are 
going to be quite energetic and can be tagged. The process pp — > ZZ — > A£ can be the 



* We discard £ — t case since such a pattern ends up with a t lepton in the final state. We concentrate 

only on the first two generations of the charged leptons. 
t Note that the r~r^ channel is not kinematically open for the LM1' scenario. 
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background for the 4£+ Et decay mode. However, a simple Et cut would eliminate events 
from the SM process pp — > ZZ — > A£. We confirmed this with our event simulation. 

At the first stage, the following basic cuts are applied to suppress the SM background. 
It is required that, whenever relevant, 

• Each isolated charged lepton (electron or muon) has a transverse momentum px(£) > 
15GeV. 

• The missing transverse energy satisfies E T > 100 GeV. 

• The leptons are constrained to be in the central barrel region of the detector by forcing 
the pseudorapidity \r)\ < 2. 

• The cone size between two charged lepton ARu is at least 0.4. Here ARu = 
(Ar) 2 + Atp 2 ) 1 ^ 2 defined in the pseudorapidity-azimuthal angle plane. 

As mention above, a missing transverse energy cut Et > 100 GeV practically gets rid of 
the SM background for the 4£+ Et signal, which is now considered background free. For 
the 4£+ Et signal, in order to get enough statistics after the cuts (as much less number of 
events pass the cuts as compared of the other two signals), we relaxed the some of the above 
cuts. We use Pt(£) > 5 GeV and AR U > 0.2 for the analysis of this signal. 

With the above cuts, the SM background is still larger than the signals 0£+ E T and 
2£+ E T . The D-Y and W + W~ dominate the ZZ cross section and they are all well above 
the signal for Et < 500 GeV. Such background domination happens in various other dis- 
tributions in most part of the region. There is no point to present these figures. Instead 
we need to find a better way to handle the background. After examining the results at 
the first stage, we decided to use E T UU1 , also known as the effective mass m eS in literature. 
This variable could be helpful in reducing the backgrounds while keeping most of the signal 
events especially if we use a suitable value for the cuts. E T um is defined as the scalar sum of 
the lepton transverse momenta and the missing transverse energy 

E s ™ = m eS = J2\PT(t)\ + $t. (19) 

l 

where the missing transverse energy E T is the sum of the total x and y components of the 
momenta in quadratures. Since it has been observed that the signal processes lead to mostly 
high m e ff (or E^ m ) distributions, a cut on m cS would substantially reduce the background. 
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Thus, as a second stage for the cuts m e g = E^ m > 750 GeV has been employed (but only 
for the 0£+ Et and 2£+ Et cases.) 

Global inclusive variables like Et and E T um are used to estimate the mass scalar of 
the parent particles produced in the hard scattering (thus estimating the scale of the new 
physics). In a recent paper by Konar et al. 23], a new global inclusive variable, called 

~l/2 

Sjnin, is proposed as an alternative. For SUSY models with R parity conservation, the decay 
chain always ends with an LSP, which is left undetected at the collider. This makes mass 
reconstruction procedure almost impossible, especially if there are more than one LSP (there 
are at least have two LSPs in the final state). Without going into extensive details of the 
signal, there is an easy way to approach guessing the scale of the new physics through the 



parameter s l J^ n . It is defined as 



23] 



y/W=I» + JV% + Visible (20) 



where E is the total calorimeter energy, P is the total visible momentum and M; nv i S ibic is the 
total mass of all invisible particles produced in each event, which is the only unknown. All 

1 /2 

the others variables can be measured at the detector. Hence s^m (^invisible) is the variable 

~l/2 

to consider. The peak of the distribution is associated with the mass threshold of the 
parent particles originated from the hard scattering. Of course, an estimation needs to be 

1 /2 

done for the total invisible mass M invisi bi e . In most of the cases the s^ in (0) gives a pretty 
good idea about the masses of the parent particles. It is shown that the method works 
better for signals with fewer invisible particles and/or more visible particles. It also works 
better with higher SUSY scales where Initial State Radiations (ISR) are less significant. We 

~ 1 /2 

include some figures for 5^(0) in the 0£+ Et as well as 4£+ Et signal. 

The events are generated at the partonic level with CalcHEP jl8| and passed to Pythia 
with the use of CalcHEP-Pythia interface for hadronization and cuts. We simulated 4 x 10 6 
events for the 0£+ Et, 2£+ Et and 4£+ Et signals. Since the relative number of events 
in each signal turns out to be proportional to the relevant branching ratio combination, 
the number of events can be simply weighed by w = a(pp — > i>(i>l) x C/N to t where £ is the 
integrated luminosity and iVtot is the total number of event generated. We set C = 100 fbT 1 , 
the ultimate goal that is expected at the LHC Even though the current reach of LHC center 
of mass energy is 7 TeV, we use 14 TeV in the numerical study, which maximizes the reach 
in the parameter space. 

19 




200 400 600 800 1000 1200 1400 1600 1800 2000 

p T /GeV 



io :i 



10- 



o 
o 



> 


o 



a 
w 



10" 



io- 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



01 + 



LM1 
LM2 
LM6 
LM1' 
LM2' 
LM6' 




200 400 600 800 1000 1200 1400 1600 1800 2000 

§Lf„(0)/GeV 




01 + Pt 



D-Y 








LM1 




LM2 




LM6 -- 




LM1' 




IA12' 




LM6' 





500 1000 1500 2000 2500 3000 3500 4000 

E^ um /GeV 




200 400 600 800 1000 1200 1400 1600 1800 2000 

8&(0)/GeV 



FIG. 3. The fi T , E^ m , and s^? (0) distributions of the 0£+ ft T signal at 14TeV with integrated 



luminosity C = 100 fb , for the three scenarios in both MSSM and the secluded U(l)' model. 



1. The Missing Energy Signal: 0£-\- $t 

The distributions of E T , and Sj£(0) are depicted in Fig. [3] for the three scenarios 

LM1', LM2' and LM6' as well as the three benchmarks for the MSSM. In general the LM6' 
scenario has the largest event pass the cuts, with similar results for the LM2', while the 
LM1' has the lowest. In fact about 80% of the events pass the cuts in LM2' and LM6' 
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but only 60% do so for the LM1'. In all three scenarios, 100% of the events pass the $t 
cut so that we are only losing 20% to 40% of them by employing the m e ff cut. This is 
because the direct LSP decay modes of the scalar neutrinos (both the left handed and right 
handed one) are either not available or suppressed for the LM1', so that the 0£+ $t signal 
would emerge from indirect decay channels through or X4 decays, with smaller branching 
ratio combinations. This can be understood from details provided earlier. The distributions 
for the secluded U(l)' model dominate the ones for the MSSM since basically the total 
production cross section in the secluded 11(1)' model is much bigger. 

The background distributions for the D-Y and ZZ processes are also included in the $t 
and E^ m graphs. The rate of success for the D-Y events passing both of the cuts are only 
about eight in 10 6 . To give an idea how effective the m e ff cut is, the success rate of events 
was about a bit more than 3% before implementing the m c g cut. The situation is even more 
drastic for the ZZ case. While the almost 100% of them passed the $t cu t, this number 
goes down to 0.3% with the m e s cut. 

We included the <v in (0) graphs in the second row of Fig. [3] to estimate the mass scale of 
the parent particles, i.e., the left handed and right handed scalar neutrinos. The graph on 
the bottom right panel is nothing but the zoom-in version of the one left handed side for 
the secluded U(l)' model. We cannot say anything about the MSSM case since the sum of 

~l/2 

the parent particle masses are varying in the 300 GeV to 600 GeV range, so that the s^ in (0) 
peak is washed out due to the m g cut at 750 GeV. Indeed, for the secluded model we should 
expect two different peaks, one for the production of the left handed scalar neutrinos and 
the other one for the right handed ones. The peak for the left-handed sneutrinos which are 
much lighter are also washed out. We will see the picture clearer for the 4£+ $t signal. 

In the right panel, we also indicated the positions of the peaks, which are correlated to 
the masses of the right handed scalar neutrinos. Of course, in reality to determine the peak 
position by fitting the data, a better job is needed . We just want to prove a point here. 
The peak position is related to the mass of vg R (since we produce them in pair) 



From the peak positions in the graph we can estimate the average right handed sneutrino 
masses ~ (530, 565, 600) GeV for the (LM1', LM2', LM6'), respectively, while the real 

average values should read (460, 563, 654) GeV from Table IIVI One source of error is not 




(21) 
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knowing the mass of the LSP (though we find out that this is not significant here since the 
LSP mass is rather light and around 100 GeV) and the other is lack of a real fitting to the 
data to pin down the location of the peaks. The estimated values are still fairly good. We 
should also note that the method works better for signals with more visible particles. 
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FIG. 4. The $t an d 2£™ m distributions of the 2£+ $t signal at 14TeV with integrated luminosity 
C = 100 fb _1 for all three scenarios in the MSSM and secluded U(l)' model. 



2. The Dilepton Signal: 2£+ # T 

We analyze the 2£+ ft? signal in a similar fashion to the 0£+ J£t on e in the previous 
subsection. The results are shown in Fig. H] and Fig. [5j The main background is from WW 
and ZZ. The D-Y does not contribute due to the transverse missing energy cut. After 
all the cuts, about 0.2% and 0.4% of the events pass for the WW and ZZ backgrounds, 
respectively. The rates were about 7% and 19%, respectively, before the m e fr cut. The 
situation for the (LM1', LM2', LM6') scenarios after all the cuts signal is (0.9%, 50%, 50%) 
survival, but 100% in each cases before the m e g cut. For the MSSM, only the LM2 gives 
2£+ signal since, for the other two scenarios, the i^e L Xi are the only final staes. The 
number of events past the cuts for the LM2 decreases to 15% from 100% after inclusion of 
the m e ff cut. 
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FIG. 5. The p T (£), M e+i - , AR e+e - and Arj e+e - distributions of the 21+ $ T signal at 14TeV with 
integrated luminosity C = 100 fb _1 for all three scenarios in both MSSM and secluded U(l)' model. 
Here, l\ represents the hardest lepton. 



We depicted the px spectra of both leptons ordered with respect to their hardness in 
Fig. |5j As expected the MSSM leptons are softer, and the distribution for ones from the 
LM2' and LM6' are very similar. The LM1' scenario is somewhere in between. In the 
invariant mass of the leptons, the LM2' and LM6' curves peak at around 60 GeV and from 
the mass spectra in Table IIV[ the mass difference 



m 7 o 

A2 



m^o is between 75 GeV to 80 GeV. 

Al 
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The 2£+ $t signal mainly goes through Xo.- For the LM1' the mass difference is 3GeV and 
not visible. The ZZ peaks at around Z boson mass as expected. 

In Fig. [5j we also include the AR,£+£- and Arje+t- distributions. It is seen that for both the 
MSSM and the secluded U(l)' model more leptons emerge with smaller separation, unlike 
the WW case the peak is at the point where the others have minimum. The background 
can be reduced further by adjusting the M e g cut value. The leptons peak when they have 
the same pseudorapidity. 

3. The Tetralepton Signal: 4£+ $t 

As we mentioned earlier, the 4£+ fix signal is practically background-free. The ZZ 
background disappears after the cut. Taking into account having relatively few 4£+ 
events, we relaxed the pt and AR cut values. It is also true that MSSM scenarios LM1, 
LM2 and LM6 do not yield a 4£+ fi T type of signal. For the LM1', only 2% of the events 
pass the cuts and among them 1% of these are 2e2/x, while the rest of the events are shared 
between 4e and 4/i. The situation is different for the LM2' and LM6'. The events which 
pass the cuts are around 68% for both cases and again half of them are the 2e2/i type and 
the rest is shared equally between 4e and 4/i. In fact, there are more 4£+ $t events in the 
LM1' scenario as compared to the other two scenarios (about 7.5% of iV tot = 4 x 10 6 for 
LMT but only 0.1% and 0.4% for LM2' and LM6', respectively). The reason is that the 
signal goes through x® which is the dominant mode for the LM1' but not for the LM2' or 
LM6'. However, the cuts reduce the LM1' events very significantly. Again the reason is the 
fact that X2 an d X? are almost degenerate for LM1', which leads to very soft leptons. 

In Fig. El fa, E^ m , and s^?(0) distributions of the 4£+ fa signal at 14TeV with 
integrated luminosity C = 100 fb _1 are shown for the three scenarios in the secluded U(l)' 

1/2 - . , . 

model. As promised, we include a Smin(0) g ra P n with the peak correlated with the V£ L u% 
production as well as the vi R v% production (no m e g cut). We can roughly tell the positions 
of the peak without doing a serious fitting. For the LM1', the first peak is around 250 GeV 
and the second one is around 950 GeV. For the LM2', they are at (470 GeV, 1070 GeV) for 
the first and the second peaks, respectively. For the LM6', the peak positions are close to 
the LM2' case, i.e., they are at (490 GeV, 1110 GeV). Then we can estimate the masses for 
LM1'/LM2'/LM6 / 

24 



ID 2 



O 

o 



> 

a 

o 



10 J 




10-1 I i i i i I i i i i I i i i i I i 

200 400 600 800 1000 1200 1400 1600 1800 2000 



p T /GeV 




200 400 600 800 1000 1200 1400 1600 1800 2000 

§Lf„(0)/GeV 




500 1000 1500 2000 2500 3000 3500 4000 

E^ um /GeV 



FIG. 6. The $ T , E^ m and s^?(0) distributions of the tetralepton (M+ $ T ) signal at 14TeV with 
integrated luminosity C = 100 fb -1 , for all three scenarios in the secluded £7(1)' model. 



m Vi ,ma I « (125 GeV, 475 GeV) / (235 GeV, 535 GeV) / (245 GeV, 555 GeV) 

V L R / est. 

while the theoretical average values obtained, including three flavors 

• (m& ,ma ) w (132 GeV, 460 GeV) / (230 GeV, 563 GeV) / (258 GeV, 654 GeV) 

V L R J theo. 

We suspect that the deviations are mainly responsible for not determining the peak position 
after fitting the data to a curve. Also, simple averaging is not quite right. One should 
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FIG. 7. The pt distribution of the 4e+ fir an d 2e2fi+ fix signals at 14TeV with integrated 
luminosity C = 100 fb _1 for all three scenarios in the secluded U{1)' model. The 4fi+ $t case is 
similar. Also the hardness of the leptons are in decreasing order. 



include a relative weight based on the relative contributions from different flavor channels. 



The pt distributions of the 4e+ $t and 2e2/i+ $t for LM6' are given in Fig. [7J The 4/x 
case is very similar to 4e. The leptons seem slightly more energetic for the 2e2/i case than 
in the other cases. LM1' and LM2' have less energetic leptons and we do not include them 
here. Fig. [8] displays two-lepton invariant mass distributions for various possibilities. As 
expected only Opposite Sign Same Flavor (OSSF) distributions have peaks at the expected 
locations since both leptons originate from the same parent unlike the other cases, Same 
Sign Same Flavor (SSSF), Same Sign Opposite Flavor (SSOF) or Opposite Sign Opposite 
Flavor (OSOF). The next figure, Fig. [9j has four-lepton invariant mass distributions for 4e 
and 2e2/i cases. The last two figures, Fig. [10] and Fig. [Til are devoted to the Ai?sB(OB) 
and A^gB(OB) distributions of the 4£+ ^ and 2e2/z+ $t signals. The subscript 'SB(OB)' 
stands for the Same Branch (Opposite Branch) and indicates where the leptons are coming 
from. We see that the distributions are very similar for At and 2e2/i. If we compare ARsb 
and Ai?QB, the former peaks at small AR while the latter peaks larger distances. For the 
pseudorapidity, even though the shape of the distributions changes, they both peak when 
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FIG. 8. Various invariant mass distributions of the Ad+ ft? signal at 14TeV with integrated 
luminosity C = 100 fb _1 for all three scenarios in the secluded U(l)' model. 



the leptons have the same pseudorapidity. 



IV. CONCLUSION 

We presented a thorough and complete analysis of the scalar neutrino production and 
decays in a U(l)' model endowed with a secluded sector. This model has several attractive 
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FIG. 9. Four lepton invariant mass distributions of 4£+ = e, /i and 2e2ii+ signals at 

14TeV with integrated luminosity C = 100 fb" 1 for all three scenarios in the secluded U(l)' model. 



features as compared to the MSSM. First, it extends the gauge symmetry to include an ex- 
tra neutral gauge boson, allowing for the presence of right handed neutrinos. Neutrinos are 
Dirac particles in this model, and masses are provided through an effective neutrino Yukawa 
coupling which is naturally suppressed by the U(l)' invariance. This model generates the /x 
term dynamically, through the VEV of a singlet scalar field. The secluded sector consists of 
three chiral superfields in addition to S, and generates correct Z'/Z mass hierarchy without 
affecting the /i parameter. Previous studies have provided extensive phenomenological anal- 
yses of this model, and notably, have provided a novel way to explain the excess positron 
flux in cosmic rays. 

The model has three right-handed scalar neutrinos, in addition to the three left-handed 
states from the SM/MSSM spectrum. Cross sections are considerably enhanced compared to 
the ones estimated in MSSM, even though for most of the parameter space studied, the signal 
is dominated by production of left handed sneutrinos, predicted to be lighter. To perform a 
through analysis, we concentrate on three MSSM benchmark parameter points, denoted by 
LMl, LM2 and LM6 and define correspondingly three U(l)' parameter points, denoted by 
LMl', LM2', LM6', specified in such that the common parameters with MSSM are identical. 
At this point, it is convenient to give a couple of remarks on the scenarios adopted here. As 
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FIG. 10. The AB,£+£- and Ar]£+£~ distributions of the 4£+ ft? signal at 14TeV with integrated 
luminosity C = 100 fb _1 for all three scenarios in the secluded U(l)' model. Here and in what 
follows, 'SB' is short-hand for 'Same Branch' and 'OB' for 'Opposite Branch'. 

pi 

we mentioned earlier, a U(l)' model with one singlet and right handed sneutrinos [8] can 
explain the excess positron flux observed by various satellite experiments. However, this 
requires a rather special mass spectrum. Indeed, it turns out that the LSP must one of the 
right handed scalar neutrinos with a mass around 100 GeV, and the next-to-LSP must also 
be a right handed scalar neutrino weighing at the TeV scale. All the other SUSY particles 
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FIG. 11. The ARg+£- and Arj£+£- distributions of the 2e2u.+ ]frj< signal at 14TeV with integrated 
luminosity C = 100 fb _1 for all three scenarios in the secluded U(i)' model. 



have to be heavier. One might ask why we did not consider such a scenario here. There 
are a couple of reasons. First of all, excess positron flux observation doesn't need to have 
an explanation coming from particle physics, only. Secondly, the signal for such a scenario 
would be mainly just missing transverse energy since all the SUSY particles other than the 
LSP are above the TeV scale so that the cross section for left-handed sneutrino production 
would be much smaller. Practically, missing energy signal with no visible particle is not 
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useful experimentally. Finally, in this work, we focused on only low-scale SUSY scenarios 
which would be discovered with the early LHC data at 14 TeV. 

After producing and decaying the sneutrinos, we identify three final-state signals: 0£+ E T: 
2£+ Et and A£+ Et and proceed to analyze them at LHC, for 14 TeV center-of-mass energy 
and with integrated luminosity C = 100 fb _1 . We compare these signals with the 0£+ Et 
and 2£+ E T signals in MSSM, and discuss the SM background (coming from Drell Yan, 
ZZ and WW production) for each. While 0£+ E T is the strongest, it has to compete with 
MSSM and suffers from considerable background suppression, while the A£+ E T signal has 
no MSSM equivalent, is practically background free, but has few events. 

We analyze the signals and suggest cuts to distinguish it from the background. In partic- 
ular E^ m = m e e, the scalar sum of the lepton transverse momenta and the missing energy 
is found to be high for the signal, thus a cut on m e s will likely reduce the background. 
Additionally a new parameter s\^ n is found to be useful for estimating the mass of parent 
particles in hard scattering. (The peak in s]^ n gives the mass threshold of left and right 
handed sneutrinos in the decay process). Using these considerations, we can estimate the 
production cross section, the products of decay and estimate the sneutrino masses. The 

~l/2 

MSSM production differs both in the number of events expected, cross section, E T nm , s^ in , 
in the 0£+ E T case; and additionally in the pr spectra of leptons (for 2£+ E T case). The 
A£+ E T case has no MSSM equivalent and little, if any, background, so the U(l)' is clear 
there; however the number of events, especially after passing detector cuts, is small. 

In conclusion, our extensive analysis shows significant enhancement of U(l)' signal over 
the MSSM signal in sneutrino production and decays, and indicates how the two models 
can be distinguished from each other and the background. This provides a distinct collider 
signal for the secluded U(l)' model at the LHC. 
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Appendix A: The Lagrangian 



In this Appendix, we present the complete Lagrangian of the U(l)' model and highlight 
the differences between this and the MSSM Lagrangian. Although parts of this formulation 
have appeared elsewhere , |25|, we include the complete model information for con- 



sistency, and to help future studies. The total Lagrangian incorporates kinetic terms and 
various interaction terms among the fields. We discuss below the distinct pieces separately. 
The kinetic terms of the Lagrangian are given by 



p Kinetic 
/ -C/(1)' 



p Kinetic 
*-MSSM 



(Al) 



where j = 1,2,3. The interactions of the gauge fields with the rest (fermions, sfermions, 
gauginos, Higgs and Higgsino fields) are contained in the piece 



pgauge pgauge ( D D w r?i \ 

L u{\)> - l mssm \9y—£>^ -> gy—t> a + 9y>Qx z h J 



(A2) 



where X runs over the fields charged under U(l)'. In (]A1|) . Z ,fJlU is the field strength tensor 
of Z'p and V u Sj = (d a + ig Y .Q' S] Z'^Sj for j = 1, 2, 3. 

The part of the U(l)' Lagrangian spanned by the F-terms is given by 



pF—term 
C/(l)' 



dW 



^mssm (a* — * h s S) h s \H u ■ Hd\ 
h„ 



- [h u Q*U* + h* s S*H*) j^S.LN 

- ^-S{L*N* (h u QU + h s SH d + ^S 1 ~LN 

Mr \ Mr 



h e H* d E*) —^-SiH u N 



Mr 

-j^-SlH*N* (h e H d E + ^-S X H U N 



h 2 



Ml 
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%\L-HJ 2 N 2 -h*S*S* 



- -tj-l* ■ h:n*ks 2 s 3 - h s s;s;^L ■ h u n 

Mr Mr 

— h s S 1 S 3 — h s S 1 S 2 (A3) 

where 0« is the scalar component of the z-the chiral superfield in the superpotential. 
The .D-term contributions to the Lagrangian are given by 

r>D—term ^ \ * r\a r^a r>D—term 



2 

,2 



'MSSM 



- 9 -f (q'qQ*Q + Q'uU*u + Q' D b*b + q' l l*l + q' e e*e 

+ Q' H H*H d + Q' Hu H:H u + Q' N N*N + Q' S S*S + Q's^S, 
+ Q' S2 S;S 2 + Q'S* 3 sX (A4) 



The soft-breaking sector of the U(l)' Lagrangian is 

C u°(iy = C-mssm^ -> 0) - m 2 s S*S - m 2 Si S* 1 S 1 - m 2 S2 S* 2 S 2 - m| 3 S*£ 3 - m 2 N N*N 
- [h s A s SH u ■ H d + ^-A V S X L ■ H U N + A- h h s S x S 2 S 3 + h.c] 

Mr 



+ -(MfrZ'Z' + h.c 

+ (mss^Si + m ss 2 SS "2 + m s 1 s 2 s i s ^ + h - c -) ( A5 ) 



where M~, is £7(1)' gaugino mass defined below in ( 1C2|) . and A s is the extra trilinear soft 
coupling. 

Finally, the part of the Lagrangian describing the fermion-sfermion-ino interactions, as 
well as the Higgs-Higgsino-Higgsino interactions, is given by 

= ^sm^ ^ 0) + iV2g Y ,[Q' Q QlZ'Q + Q' v u\Z'ur 

+ Q' D d R Z'd R + Q' L L^Z'L + Q' E £ R Z'£ R + Q H E\Z'E d 

+ Q' Hu HlZ'H u + $ 8 &Z'S + Q' Sj S]Z'S j + Q' N u R Z'u R + /i.e.] 

+ [h s SH u ■ H d + h s SH u ■ H d + /i s £# u • H d + kc.]. (A6) 

All parts of the the U(l)' model Lagrangian listed above are described in the current 
basis. Eventually, the fields must be transformed into the physical basis where each field 
obtains a definite mass. The neutral gauginos and Higgsinos form the neutralino sector 



whose physical states are expressed as in (]C1|) . after diagonalizing the mass matrix flC2j) 
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Unlike the neutralino sector, the structure of the chargino sector is essentially the same as 
in the MSSM with the replacement /i — > h s v s /\/2. A detailed analysis of the Higgs and 
chargino sectors of the U(l)' model has been given in 

In the gauge boson sector, spontaneous breakdown of the product group SU(2) l®U(1)y® 
U(l)' via the Higgs VEVs 

(Hu) = ^= \ ° I (//,/>■ ~= I I ($)■■ -7=- <*',-) (AT) 




generates one massless state (the photon) and two massive states (the Z, Z' bosons) via 
orthonormal combinations of W^, B' and B^ gauge bosons. The and W 2 linearly 
combine to give W^, as the only charged vector bosons in the model. In contrast to the 
MSSM, the Z boson is not a physical state by itself since it mixes with the Z' boson. This 
mass mixing arises from the fact that the Higgs doublets H U ( i are charged under each factor 



of SU{2)i J £g> U(l)y <8> U(l)', and the associated mass-squared matrix is given by 

Mf A 2 
A 2 Mf, 
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, Ml A 2 , 

ML* = ! „ > (A8) 



in the Z') basis. Its entries are 

Mf = \g\ {yl + v 2 ) 



mi = & Q' H y u + Q' H y d + q'/v 2 s + y: & 



V 



A 2 = \G z g Y , (Q' H y u - Q' H v 2 d ) , (A9) 

where G 2 Z = g^ + Qy- The physical neutral vector bosons, Zip, are obtained by diagonalizing 
M 2 z ^ z r. 



- sin 8 Z -z> cos 6 z-z 




(A10) 



where 

fl„ „, = _ <(J , i<( ,, , 

2 V M 2 - Ml 



1 / 2A 2 
9z-z' = — -arctan — = — j ) , (All) 



is their mass mixing angle, and 



M i( 2 ) = \ 



M\, + M 2 Z - (+) V(M 2 , - M 2 ) 2 + 4A 4 



(A12) 
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are their masses-squared . The collider searches at LEP and Tevatron plus various indirect 



observations require Z—Z' mixing angle Qz-z 1 to 
model dependence coming from the Z' couplings 



3e at most a few 10 3 with an unavoidable 
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31] . This bound requires either Mz 2 to 



be large enough (well in the TeV range) or A 2 to be sufficiently suppressed by the vacuum 
configuration, that is, tan 2 (3 = v 2 Jv 2 d ~ Q'h IQ'h u - Which of these options is realized 
depends on the U(l)' charge assignments and the soft-breaking masses in the Higgs sector 
( see [3| for a variant for reducing the Z—Z' mixing). 



Appendix B: The Scalar Fermions 

Given rather tight FCNC bounds, we neglect all the inter-generational mixings, and con- 
sider only intra-generational left-right mixings, though these turn out to be totally negligible 
for the sfermions in the first and second generations. The 2x2 scalar fermion mixing matrix 
can be written as 



/ 



M 2 ~ v, , 

* 1 fa J 1 io.,b 
■ILL J LR 



M% b M 2 ~ 

J LR J R 



a 7^ b = u,d . 



(Bl) 



where 



M) l + \h%vlK 2 s + i 



9yYf ? 



+ ^9Y>Q'f?{Q'H u Vu + Q'H d V d + Q's V sP 



( + )| K-vi) 

I „.2 



Ms 

JR 



1 il 2 2 | 1 



+-g Y >Q'f«{Q'H u v u + Q r H d v d + Q's v sP 



M% >l3 = (M%, fi ) 

Jlr Jrl 



h fa K s 



where k. 



V2M R 



and £ s 



RL 

h s v sn v 3 „v^ 



(±2A* fa v a + V2h s V0V s + 2V2£s 



(B2) 

(B3) 
(B4) 



2^2 v "^ /Q 

. ua 2^3 f or sri eutrinos and k s = 1 and £ s = for the others. 



Here M~ are the soft mass-squared of the sfermions, v Ui ^ tStSljS2iS3 are the VEVs of the Higgs 

jL,Ii 

fields, Yfa(T 3L ) is the U(l)y (SU(2) L ) quantum number, Q'* a is the U{1)' charge, and Afa 
are the trilinear couplings. The mixing matrix can be diagonalized, in general, by a unitary 
matrix such that V*' ■ Ai 2 ~ ■ V*" = Diag(M~ a , M~ a )H| The rotation matrix r*° can be 

* We note that unlike mixings in other sectors, is defined differently, that is, (/£ R )i — t{j /?, where 
fj represent the mass eigenstates. 
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written for quarks and charged leptons in the 2x2 {/£, /^} basis as 



fa _ 1 cos 9j a - sin 9j a 
sm6j a cos9j a 



(B5) 



where 9j a 



— arctan2(— 2M. \ a ,-Mr a — M 2 - ) and arctan2(y, x) is defined as 

2 JLR J RR J LL 



arctan2(y, x) 



4> sign(y), x > 

f sign(y), x = 

(tt - 0) sign(y), x<0 



(B6) 



with y being non-zero, and <fi taken in the first quadrant such that tan0 = \y/x\. 

For the sfermions in the first and second generations, the left-right mixings are exceedingly 
small as they are proportional to the corresponding fermion mass. Therefore, the sfermion 
mass matrix (1B2|) is automatically diagonal. However, one has to remember that the sfermion 

masses, for fixed values of ms , are different in the MSSM than in the U(lY models due 

Jl,r 

to the additional D-term contribution in the latter. 



Appendix C: Gauge and Higgs Fermions 

Although the U(l)' model possesses no new charged Higgsinos and gauginos it possesses 
five new fermion fields in the neutral sector: the U(l)' gauge fermion Z' and four singlinos 
S, Si, S2, S3. In total, there are 9 neutralino states %° (i — 1, . . . , 9) 

5ff = XXG a , (CI) 



where the mixing matrix N® a connects the gauge-basis neutral fermion states G a G | -B, W 3 , 
H%, H®, S, Z', Si, S2, S3 1 to the physical neutralinos Xi- The neutralino masses M^o and 
the mixing matrix iV? are determined via the diagonalization condition N°A4N° T = Diag 



M^o , . . . , M^o I for the neutral fermion mass matrix 
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(C2) 



where certain entries are generated by the soft-breaking sector while others follow from the 
SU(3) C <g> SU(2) L ® C/(l)y <g> C/(l)' breaking. The I7(l)y gaugino mass M y , the 5C/(2) Z 
gaugino mass M^, and the U(l)' gaugino mass 



My, tan y „ , 9 

> -2 — M yy , + My tan x 



cos 2 % 



cosx 



as well as the mixing mass parameter between £/(l)y and U(l)' gauginos 

M vv , 

M y ^, = -Mptanx, 



(C3) 



(C4) 



all follow from the soft-breaking sector. Through the mixing of the gauge bosons, M z , and 
Myg, exhibit an explicit dependence on the masses of the U(l)y and U(l)' gauginos, and 
their mass mixing. M yy , is the soft-breaking mass that mixes the U(l)y and U(l)' gauginos. 
In the numerical analysis, we set the mixing mass parameter M YZ , = since we neglect the 
kinetic mixing (tan x — > 0) thus M YY , — > 0. For convenience we also define Ry = My, /My. 

The remaining entries in flC2j) are generated by the soft-breaking masses in the Higgs 
sector via the SU(3) C ® SU(2)l <g> U(l)y <8> C/(l)' breaking. Their explicit expressions are 
given by 

My g d = M z sin 6 W cos /3 , My = M z sin 6 W sin /3 , 



M^? g = M z cos 6^ cos ft , ^ = M z cos 6^ sin ft , 

v'h u = 9y'Q' Hu v u , /% = gvQ'sVs , a4 = gyQ's^si , 
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(C5) 



where gy is the coupling constant of U(l)'. For numerical analysis we choose the standard 
GUT value for it gy = J^gtax^Ow 

Appendix D: The Compositions of the Neutralinos 

In this Appendix we give the Bino, Wino, Higgsino and Singlino compositions of the 
physical neutralinos = 1,2,..., 9 for the three scenarios LM1', LM2' and LM6'. They 
are listed in Table IVII 
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LMl' 


Xi 


X2 


X°3 


xl 


xl 


xl 


X? 


x^ 


x|j 


B 


-0.988 


0.046 


0.077 


0.043 


-0.056 


0.095 


-0.002 


0.0 


0.0 


W 3 


0.037 


-0.058 


0.955 


-0.122 


0.086 


-0.245 


0.006 


0.0 


0.0 


H° d 


-0.126 


0.051 


-0.245 


-0.321 


0.692 


-0.581 


0.014 


0.0 


0.003 


H° u 


0.031 


-0.205 


0.115 


0.226 


0.698 


0.633 


-0.033 


0.002 


-0.040 


S 


0.057 


0.910 


0.086 


0.362 


0.140 


0.0 


0.048 


-0.004 


0.087 


Z' 


-0.010 


-0.180 


-0.010 


0.212 


0.022 


-0.092 


0.065 


0.013 


0.953 


Si 


-0.012 


-0.148 


-0.013 


0.492 


0.012 


-0.276 


-0.576 


0.555 


-0.133 


s 2 


-0.006 


-0.089 


-0.006 


0.155 


0.008 


-0.060 


0.778 


0.586 


-0.120 


§3 


0.018 


0.241 


0.020 


-0.621 


-0.023 


0.320 


-0.232 


0.589 


0.223 


LM2' 


B 


0.048 


-0.994 


0.032 


-0.009 


-0.044 


0.076 


-0.005 


0.0 


0.0 


W 3 


-0.054 


0.0101 


0.974 


0.029 


0.067 


-0.204 


0.014 


0.0 


0.0 


H° d 


0.032 


-0.088 


-0.191 


0.135 


0.699 


-0.666 


0.043 


0.0 


0.001 


Hi 


-0.224 


0.013 


0.083 


-0.062 


0.692 


0.673 


-0.057 


0.002 


-0.033 


S 


0.942 


0.054 


0.075 


-0.238 


0.152 


0.130 


0.028 


-0.003 


0.069 


Z' 


-0.119 


-0.006 


-0.007 


-0.190 


0.015 


-0.019 


0.047 


0.011 


0.972 


Si 


-0.091 


-0.007 


-0.010 


-0.550 


0.006 


-0.143 


-0.585 


0.561 


-0.100 


s 2 


-0.069 


-0.005 


-0.006 


-0.212 


0.005 


0.013 


0.778 


0.578 


-0.094 


S3 


0.165 


0.012 


0.016 


0.732 


-0.013 


0.123 


-0.206 


0.591 


0.169 



LM6' 



B 


-0.036 


-0.995 


-0.034 


-0.003 


0.035 


-0.072 


0.006 


0.0 


0.0 


W 3 


0.041 


0.015 


-0.978 


0.014 


-0.054 


0.194 


-0.017 


0.0 


0.0 


H° d 


-0.015 


-0.080 


0.174 


0.093 


-0.700 


0.677 


-0.060 


0.0 


-0.002 


H u 


0.196 


0.021 


-0.089 


-0.015 


-0.695 


-0.680 


0.073 


0.001 


0.031 


S 


-0.936 


0.040 


-0.057 


-0.278 


-0.143 


-0.114 


-0.028 


-0.004 


-0.073 


z' 


0.138 


-0.005 


-0.006 


-0.202 


-0.012 


0.007 


-0.046 


0.012 


-0.967 


Si 


0.116 


-0.007 


0.010 


-0.541 


-0.006 


0.120 


0.593 


0.560 


0.109 


S2 


0.082 


-0.004 


0.005 


-0.215 


-0.004 


-0.044 


-0.772 


0.581 


0.101 


S3 


-0.203 


0.011 


-0.015 


0.729 


0.012 


-0.070 


0.197 


0.589 


-0.185 



TABLE VI. The Bino, Wino, Higgsino and Singlino composition of the neutralinos x®, i = 1, 2, 9 
for the scenarios LM1',LM2' and LM6'. 
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